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Abstract 

The compound A — {BETS)2FeCli provides an effective demonstration of the interaction of 
TT conduction electron and d-electron localized moment systems in molecular crystalline materi- 
als where antiferromagnetic insulating and magnetic field induced superconducting states can be 
realized. The metal-insulator transition has been thought to be cooperative, involving both the 
itinerant vr- electron and localized d-electron spins where antiferromagnetic order appears in both 
systems simultaneously. However, recent specific heat data has indicated otherwise [Akiba et al, J. 
Phys. Soc. Japan 78,033601(2009)]: although the vr-electron system orders antiferromagnetically 
and produces a metal-insulator transition, a "mysterious" paramagnetic d-electron state remains. 
We report ^'^ Fe Mossbauer measurements that support the paramagnetic model, provided the d- 
electron spins remain in a fast relaxation state below the transition. From the measured hyperfine 
fields, we also determine the temperature dependence of the n — d electron exchange field. 

PACS numbers: 74.70.-b, 74.25.Bt, 74.70.Ad 
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A — {BETS)2FeCl4^ (BETS = bisethylenedithio-tetraselenafulvalene) is one of the most 
thoroughly studied molecular conductors in the last few years due to its unique proper- 
ties derived from the interaction between conducting vr-electrons in the BETS donor lay- 
ers and localized d-electrons in FeC/4 anions with S=5/2 spins. ^- The crystal structure of 
A — {BETS)2FeCli as shown in Fig. 1 consists in stacks of BETS donors along a and tightly 
packed in layers with a 2D network of S...S contacts parallel to (a, c) , alternating along h 
with layers of FeCl^ anions. At high temperatures this compound is a quasi 2D metal due 
to delocalization of vr-electrons in the layered network of partially oxidized donors, and its 
magnetic susceptibility is dominated by the paramagnetic S=5/2 FeCU spins. At 8.3 K 
this compound undergoes a transition towards an antiferromagnetic (AF) insulating ground 
state (see inset of Fig. 3 below). Since the isomorphous compound with diamagnetic GaC/4 
anions remains metallic, becoming superconductor at 6 K, the metal-insulator transition in 
A — {BETS)2FeCl/i has been thought to be driven by the ordering of the Fe spins. 

However until now no direct microscopic measurements directly probing the role of the 
anions have been published and recent specific heat measurements by Akiba et al. have cast 
some doubts on the role of the S=5/2 FeCU spins in the transition.- These authors have 
suggested that during the transition, while the vr spins order antiferromagnetically, the Fe 
spins remain paramagnetic below 8.3K. According to their model an effective field i^Tr-d ~ 4 
T caused by the ordering of the vr spin system is switched on at the Fe sites at approximately 
8.3 K but the Fe^^ cations remain paramagnetic with the 3d energy levels described by a 
Zeeman splitting. The latter gives rise to a Schottky 6- level term in the specific heat. 

In this Communication we describe the results of a ^''Fe Mossbauer spectroscopy study 
to examine the role of the Fe S=5/2 spins in the transition. The single crystals of 
A — {B ET S)2F eCli used in this work were grown using standard electrochemical meth- 
ods from 99% ^'^Fe enriched TEAFeCli. Mossbauer spectra were collected with the 
absorber within a liquid-He bath cryostat, in transmission mode using a conventional 
constant-acceleration spectrometer and a 25 mCi ^^Co source in a Rh matrix. The ab- 
sorber was prepared by randomly placing between two perspex plates approximately 4 mg 
of A — {B ET S)2F eCli single crystals 99% enriched in ^'^ Fe. Electrical transport measure- 
ments, which verified the metal-insulator transition in the ^''Fe enriched samples at T/v = 
8.3 K (inset of Fig. 3), were carried out using a standard 4-terminal resistance configuration 
on crystals that had been used in the Mossbauer measurements. 
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FIG. 1: Crystal structure of \ — {BETS)2FeCli projected in the b-c plane(after Ref.-). The charge 
transfer of one electron between each two BETS donors and one FeCl^ anion leaves a delocalized 
spin 1/2 TT electron on the donor stacks, and a localized d-electron at the anion site. The schematic 
shows the two possible spin configurations below T/v with fully antiferromagnetic tt and d order 
(solid arrows) or with a paramagnetic d state (dashed arrows). 

Mossbauer spectra obtained at different temperatures in the range 9 to 1.5 K are shown in 
Fig. 2, and the fitting parameters are listed in Table 1. The spectra were fitted to Lorentzian 
lines using a non-linear least-squares method.™ The relative areas and line widths of both 
peaks in a quadrupole doublet and the peak pairs 1-6, 2-5 and 3-4 in a magnetic sextet 
were constrained to remain equal during the refinement procedure. Isomer shifts (IS) are 
given relative to a — Fe at room temperature. Above 8.3 K the spectra present a single 
line typical of paramagnetic Fe^~^ in a site with a low quadrupole distortion (high symmetry 
environment). Below 8.3 K sextet splittings are observed. The spectrum obtained at 8 K 
still shows a small fraction (ca. 13%) of "paramagnetic" atoms, which can be ascribed to 
a hysteresis or the slow dynamics of the transition. Below this temperature, down to 3.2 
K, two magnetic splittings are clearly seen. They have identical isomer shifts, but slightly 
different hyperfine fields, B^f . This clearly indicates that there are two different magnetic 
environments for the Fe sites. The uncertainty in the relative areas of the sextets is relatively 
large, but in a first approximation they present a temperature independent 1:2 ratio. Below 
3 K, these two sextets appear to merge and only one sextet is observed in the range 2.8-1.5 
K, suggesting a transition occurring at 3.0 ± O.IK. The hyperfine fields Bhf , shown in Fig. 
3, present significant temperature dependence until circa 3 K. Below this temperature they 
seem to saturate approaching 44 T at 1.5 K, very close to the value of B^f =45 T observed 
for FeCh at 4.2 K.^ 

The observed change of the shape of the Mossbauer spectra in a narrow temperature 
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FIG. 2: Fe Mossbauer spectra of A — {BETS)2FeCU for different temperatures. The indi- 
vidual solid lines are components of the spectrum coming from the two sextets and quadrupole 
components. The solid line through the data (dots) is the sum of the individual contributions. 

range around 8.3 K is indicative of a magnetic ordering process. Clearly, above the M-I 
transition the Fe^'^ spins are oscillating with a relaxation frequency u^i > 10^^ Hz, typical 
of paramagnetic Fe^^. In this situation the magnetic field observed at the Fe nuclei is 
averaged to zero.-"^ In a first approximation the presence of sextets with sharp peaks as 
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FIG. 3: Temperature dependence of the Fe hyperfine fields, B^f in A — {BETS)2FeCli . The 
lines represent the fast relaxation model for both fixed and temperature dependent Ht^_ci exchange 
fields (see text). Inset: temperature dependent resistance of a single crystal of ^"^Fe enriched 
A - {BETS)2FeCU showing the M-I transition at 8.3 K. 

those observed below the M-I transition suggests that the Fe nuclei are feeling a static 
magnetic hyperfine field within the observation time scale of the Mossbauer effect. This 
would imply that the relaxation frequency of the Fe^^ spins is now ujr < 10^ Hz. 

In the proposed picture of Akiba et al.- based on specific heat data of \ — {BETS)2FeCU , 
the Fe S=5/2 spins remain paramagnetic below the transition at 8.3 K, although subject 
to an internal effective field H-j^^d ~ 4 T switched on at the Fe sites by the ordering of 
the TT spin system. In the absence of -ffjr-d; if the Fe atoms remain in the paramagnetic 
state in a fast relaxation regime, the magnetic field in the nuclei averages to zero and only 
a two-line pattern should appear corresponding to the electric quadrupole interaction (in 
this case a single absorption peak is observed due to the very low quadrupole splitting), 
evident at 9 K. If there was a drastic slowing down of the electronic relaxation, but in the 
absence of magnetic ordering, two different sextets should be observed, corresponding to the 
different Ms states 5/2, and 3/2, as it has been observed in diluted ^^Fe doped compounds, 
such as in Fe-doped Al203,—LiAl50^ or in proteins such as the transferrins.™ The M5=l/2 
state should result in a complicated 11-line pattern since it induces non-diagonal terms in 
the Hamiltonian of the hyperfine interactions.-^ Although its presence was reported in the 
transferrin case^ it usually is not observed due to enhanced relaxation of the Ms = ± 1/2 
electronic doublet .- 
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Unlike dilute systems, the concentration of Fe^^ in A — {BETS)2FeCU is not so low as 
in the above examples and therefore due to spin-spin interactions a slow relaxation regime is 
not expected to occur. It could be argued that in A — {BETS)2FeCli the internal effective 
field of the donors at the Fe sites, estimated as H^^^d ~ 4 T, could freeze the spin flipping 
and bring the system to a slow relaxation regime. However there is no evidence for the Ms 
= 3/2 state. The sextet with the smaller magnetic splitting observed between 8 and 3K 
cannot correspond to the Ms = 3/2 state because its Bhf value is > 90% of the Bhj of 
the larger magnetic splitting (below 4.1 K where thermal excitations are less important). 
The saturated Bhf values associated with Ms =3/2 and 5/2 electronic states should be 
proportional to Ms-^ Furthermore in such case the relative intensities of both magnetic 
splittings should follow the evolution of their statistical thermal population, while they 
remain approximately constant. 

The origin of the transition observed at approximately 3 K, as the merging of the two 
sextets with only one hyperfine field is not entirely clear, but may be associated with a 
change of the magnetic wave vector in the AF state. Although not as dramatic as in the 
Mossbauer data, changes below T^r have been seen in other independent studies in the 3 
K region. Matsui and co-workers^i have investigated the microwave cavity response with 
A — {BETS)2FeCli , and below have found highly dispersive modes attributed to charge 
degrees of freedom. However, for Hac \\ H || a*, a peak in the cavity dissipation (Ar/2/o) 
appears at 3 K (H = 0) which, due to the unfavorable direction of the eddy currents for 
Hac II a*, the authors attribute to a dynamic response due to spin degrees of freedom. This 
peak has a complicated dependence on H near the 1.2 T spin-flop transition. Likewise, Rutel 
et al.— have observed anomalies in the microwave cavity response below 4 K for H || c— 
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Table I. Computed parameters from the Mossbauer spectra of A — {BETS)2FeCU taken 
at different temperatures. 
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IS 


QS, e 


Bhf 


r 


I 


(K) 


(mm/s) 


(mm/s) 


(T) 


(mm/s) 




200 


0.27 


0.21 


- 


0.49 


100% 


9.0 


0.34 


0.20 


- 


0.54 


100% 


8.0 


0.34 


0.20 


- 


0.57 


13±2% 




0.33 


0.10 


17.3 


0.33 


26±4% 




0.33 


0.07 


22.9 


0.40 


61±4% 


7.2 


0.34 


0.11 


20.1 


0.41 


32±3% 




0.34 


0.08 


26.0 


0.39 


68±3% 


6.1 


0.34 


0.10 


23.9 


0.47 


32±2% 




0.34 


0.10 


29.6 


0.52 


68±2% 


4.1 


0.34 


0.08 


33.5 


0.43 


36±1% 




0.34 


0.09 


37.6 


0.42 


64±1% 


3.7 


0.34 


0.07 


35.3 


0.29 


43±3% 




0.34 


0.08 


38.8 


0.33 


57±3% 


3.25 


0.33 


0.08 


37.6 


0.32 


38±1% 




0.33 


0.08 


40.6 


0.36 


62±1% 


2.8 


0.33 


0.09 


41.2 


0.45 


100% 


2.4 


0.33 


0.07 


42.5 


0.38 


100% 


1.8 


0.34 


0.09 


43.7 


0.61 


100% 


1.5 


0.34 


0.09 


44.0 


0.44 


100% 



IS, isomer shift relative to metallic Fe at 298 K. 

QS, quadrupole splitting. 

e — {e^VzzQ / ^){^cos^9 — 1), quadrupole shift. 

Bhf, magnetic hyperfine field. 

r, half-width of the doublet peaks. 

I, relative area. 

Estimated errors: < 0.002 mm/s for IS, QS, e, T; < 0.2 T for Bhf 

It is fiowever difficult for a Mossbauer probe to discriminate between tfie onset of magnetic 
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order, spin-glass, behavior, or a particular case of a paramagnetic "fast relaxation" behavior 
where the Fe spins are Zeeman split by an applied field and the population of M5 states 
are different. The present results may therefore be consistent with the Fe atoms remaining 
paramagnetic in the low temperature state below 8.3 K. Assuming a fast relaxation model 
below the transition, we may compute the hyperfine field based on the Fe cations, inducing 
a Zeeman splitting and a Boltzmann distribution of the 6 states: 

Bhf{T) = Y.[B{Ms)exp{XMs)]/Y.^xp{XMs) (1) 

Ms Ms 

Here Xm, = -gfiBMsH^-d/kBT. We take B{M,) ~ ±45 T for M,=t5/2, ±27 T for 
Ms=^3/2, and ±9 T for Ms==i=l/2 (in proportion to Mg). In Fig. 3 the computed temper- 
ature dependence of B^f is shown for both a constant H^^^^ = 4 T, and for a temperature 
dependent Hj^^d, increasing from 2.45 T at 8 K to 4.2 T at 1.5 K. We find that the best fit 
to the data implies that -ffjr-d is temperature dependent. Although this modifies the tem- 
perature dependence of the Schottky specfic heat described by Akiba et al.'^, the differences 
are not significant. 

The temperature dependence of the exchange field may be described by a spin-wave 
behavior—!^ where, for anti-ferromagnetic dispersion in the spin 1/2 vr-electron system uj = 
J/h\ka\, H^^d{T)= if^_d(0)(l - AT'^). We find that for A ~ 7 x IQ~^K-^, Eq. 1 provides 
a reasonable description of the temperature dependence of -B^/, shown as the solid line in 
Fig. 3. From A we estimate the spin-wave exchange energy to be J ~ 5.6/^, comparable 
to T/v. 

The suggestion and evidence that the magnetic order appears in the vr-electron system, 
but not in the d-electron system, seems unusual^. However, estimates do show that the 
mean-field exchange interaction of the vr-electron system is the largest: J-n--Ki Jn-d, and 
Jd-d are 448, 14.6, and 0.64 K respectively^. Nevertheless, the d-electron spins must play 
a central role in the formation of the magnetic ground state. A temporal probe of the 
spin dynamics of the d-electron system below the M-I transition, as well as magnetic field 
dependent specific heat and Mossbauer experiments, would be useful to further explore the 
nature of the magnetic order associated with this very unusual antiferromagnetic metal- 
insulator transition. 



8 



Acknowledgments 



This work was supported by FCT(Portugal) through project PTDC/QUI/64967/2006 , 
and in part by NSF DMR-0602859 (JSB). The National High Magnetic Field Laboratory is 
supported by NSF DMR-0654118, by the State of Florida, and the DOE. The collaboration 
between Portuguese and American team members was supported by FLAD. The authors 
are indebted to Pedro Schlottmann and Kun Yang for valuable discussions. 



1 H. Kobayashi, H. Cui and A. Kobayashi, Chem. Rev. 104, 5265 (2004). 

^ H. Kobayashi, H. Tomita, T. Naito, A. Kobayashi, F. Sakai, T. Watanabe, P. Cassoux, J. Am. 

Chem. Soc. , 118, 368(1996). 
^ H. Akiba, S. Nakano, Y. Nishio, K. Kajita, B. Zhou, A. Kobayashi and H. Kobayashi, J. Phys. 

Soc. Japan 78, 033601 (2009). 
^ J. V. Rodrigues, I. C. Santos, V. Gama, R. T. Henriques, J. C. Waerenborgh, M. T. Duarte 

and M. Almeida, J. Chem. Soc. Dalton Trans. 18, 2655 (1994). 
^ N. N. Greenwood and T. C. Gibb, Mossbauer Spectroscopy (Chapman and Hall, Ltd. , London, 

1971). 

^ H. H. Wickman and G.K.Wertheim, in Fundamentals of Gamma-Resonance spectroscopy, edited 
by V. L Goldanskii and R. H. Herber (Academic Press, New York, 1968), p. 554. 
H. H. Wickman and G. K. Wertheim, Phys. Rev. 148, 211 (1966). 

^ P. J. Viccaro, F. d. S. Barros and W. T. Oosterhuis, Phys. Rev. B 5, 4257 (1972). 

^ S. A. Kretchmar, M. Teixeira, B. H. Huynh and K. N. Raymond, Biol. Metals 1, 26 (1988). 

° H. D. Pfannes, R. Magalhes-Paniago and A. Putzka, Hyperf. Literac. 70, 1133 (1992). 

^ H. Matsui, H. Tsuchiya, T. Suzuki, E. Negishi and N. Toyota, Phys. Rev. B 68, 155105 (2003). 

^ L B. Rutel, S. Okubo, J. S. Brooks, E. Jobiliong, H. Kobayashi, A. Kobayashi and H. Tanaka, 
Phys. Rev. B 68, 144435 (2003). 

3 C. Kittel, Introduction to Sohd State Physics, 3^'^ Ed. (Wiley New York, 1967). 

^ Here we assume that Hpi^^ at the Fe site follows the fractional change of staggard magnetiza- 
tion in the vr-electron system. Following Refji^ for linear antiferromagnetic dispersion, the 
prefactor is A = 2A0411 {In^ / SQ){kB /2JSf where S, J, and Q are the spin, exchange energy. 



9 



and Fe sites per unit cell respectively, and the pre-factor is the value of the definite integral. 

S. Uji, J. Phys. Soc. Jpn. 78, 033601 (2009). 

H. Mori and M. Katsuhara, J. Phys.Soc. Japan 71, 826 (2002). 



10 



